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A limited mobility nonequilibrium solid-on-solid dynamical model for kinetic surface growth is 
introduced as a simple description for the morphological evolution of a growing interface under 
random vapor deposition and surface diffusion bias conditions. Large scale stochastic Monte Carlo 
simulations using a local coordination dependent instantaneous relaxation of the deposited atoms 
produce complex surface morphologies whose dynamical evolution is not consistent with any of 
the existing continuum dynamical surface growth equations. Critical exponents for coarsening and 
roughening dynamics of the morphological evolution are quantitatively calculated using large scale 
simulations. 
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An atom moving on a free surface is known to en- 
counter an additional potential barrier, often called a 
surface diffusion bias JlJ (or an Ehrlich |g| - Schwoebel ||] 
barrier), as it approaches a step from the upper terrace 
- there is no such extra barrier for an atom approaching 
the step from the lower terrace (the surface step sepa- 
rates the upper and the lower terraces). Since this diffu- 
sion bias makes it preferentially more likely for an atom 
to attach itself to the upper terrace than the lower one, 
it leads to mound (or pyramid) - type structures on the 
surface under growth conditions as deposited atoms are 
probabilistically less able to come down from upper to 
lower terraces. This dynamical growth behavior under a 
surface diffusion bias is sometimes called an "instability" 
because a flat two dimensional surface growing under a 
strong surface diffusion bias is unstable toward three di- 
mensional mound/pyramid formation. There has been a 
great deal of recent interest |o| in the morphologi- 
cal evolution of growing interfaces under nonequilibrium 
growth conditions in the presence of such surface diffu- 
sion bias. In this Letter we introduce, what we believe to 
be, the minimal nonequilibrium atomistic growth model 
for ideal molecular beam epitaxial - type random vapor 
deposition growth under a surface diffusion bias. Exten- 
sive stochastic simulation results presented in this paper 
establish the morphological evolution of a surface grow- 
ing under diffusion bias conditions to be rather complex. 
Various critical growth exponents pl] , p^| , which asymp- 
totically describe the large-scale dynamical evolution of 
the kinetically growing surface in our minimal discrete 
growth model, are inconsistent with all the proposed con- 
tinuum theories for nonequilibrium surface growth under 
diffusion bias conditions. Our results lead to the con- 
clusion that a continuum description for nonequilibrium 
growth under a surface diffusion bias does not exist and 
may require a theoretical formulation which is substan- 
tially different from the ones currently existing in the 
literature. 

In Fig. 1(a) we schematically describe our solid-on- 
solid (SOS) nonequilibrium growth model : (1) Atoms 
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FIG. 1. (a) Schematic configuration defining growth rules 
in 1+1 dimensions, (b) Morphologies (Pu = l;P L = 0.5) for 
L = 100 x 100 at 10 3 ML and (c) 10 6 ML. 
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are deposited randomly (with an average rate of 1 
layer/unit time, which defines the unit of time in the 
growth problem - the length unit is the lattice spacing 
taken to be the same along the substrate plane and the 
growth direction) and sequentially on the surface start- 
ing with a flat substrate; (2) a deposited atom is incorpo- 
rated instantaneously if it has at least one lateral nearest 
- neighbor atom (i.e. if it has a coordination of 2 or 
more since there must always be an atom underneath to 
satisfy the SOS constraint); (3) singly coordinated de- 
posited atoms (i.e. the ones without any lateral neigh- 
bors) could instantaneously relax to a neighboring site 
within a diffusion length of I provided the neighboring 
site of incorporation has a higher coordination than the 
original deposition site; (4) the instantaneous relaxation 
process is constrained by two probabilities Pl and P\j 
(with < Pl,Pjj < 1) where Pl(u) is the probability for 
the atom to attach itself to the lower(upper) terrace after 
relaxation (note that a "terrace" here could be just one 
other atom). The surface diffusion bias is implemented 
in our model by simply taking Pjj > Pl , making it more 
likely for atoms to attach to the upper terrace. Under 
the surface diffusion bias, therefore, an atom deposited 
at the top of a step edge feels a barrier (whose strength 
is controlled by Pu — Pl) in coming down compared with 
an atom at the lower terrace attaching itself to the step. 
Our surface diffusion bias model is well-defined for any 
value of the diffusion length I including the most com- 
monly studied situation of nearest - neighbor relaxation 
(/ = 1). We have carried out extensive simulations both 
in 1 + 1 and 2+1 dimensions (d) varying Pl, Pjj as well as 
I, also including in our simulations the inverse situation 
(the so-called 'negative' bias condition) with Pl > Pjj so 
that deposited atoms preferentially come down attach- 
ing themselves to lower steps producing in the process a 
smooth growth morphology. Because of lack of space we 
do not present our negative bias results here except to 
note that the smooth dynamical growth morphology un- 
der our negative bias model obeys the generic Edwards - 
Wilkinson universality jl 21 ,22]]. 

Before presenting our numerical results we point out 
two important features of our growth model : (1) For 
Pl = Pjj = 1 our model reduces to the one introduced in 
ref. 23 (and studied extensively pH-|27|| in the literature) 
as a minimal model for molecular beam epitaxy in the 
absence of any diffusion bias; (2) we find, in complete 
agreement with earlier findings p3|,^4| in the absence of 
diffusion bias, that the diffusion length I is an irrelevant 
variable which does not affect any of our calculated crit- 
ical growth exponents (but does affect finite size correc- 
tions - increasing I requires a concomitant increase in the 
system size to reduce finite size effects). We, therefore, 
mostly present our / = 1 simulation results here empha- 
sizing that our critical exponents are independent of I 
provided finite size effects are appropriately accounted 
for. Our calculated exponents are also independent of 
the precise values of Pjj and Pl as found in ref. 23 and 
24 for the P\j = Pl case. 
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FIG. 2. (a) The mound evolution in 1+1 dimensions 
showing a section of 500 middle lattice sites from a sub- 
strate size of L = 10000 at 10 4 ML, 10 5 ML, and 10 6 ML 
(Pu = 1;P L = 0.5). (b,c) The d=2+l growth morphologies 
on a 100 x 100 substrate (b) 10 3 ML, and (c) 10 6 ML. The 
darker (lighter) shades represent lower (higher) points on the 
surface. 



In Fig.l we show representative d=2+l simulated 
growth morphologies for our positive bias model. The 
diffusion bias produces mounded structures which are vi- 
sually statistically scale invariant only on length scales 
much larger (or much smaller) than the typical mound 
size. In Fig. 2 we show morphological evolutions for 
d=2+l and 1+1 in order to visualize the coarsen- 
ing/steepening dynamics of the mound structures. We 
mention that in producing our final results we utilize a 
noise reduction [ p8| technique which accepts only a frac- 
tion of the attempted kinetic events, and in the process 
produces smoother results (reducing noise effects) with- 
out in any way affecting the critical exponents. The cor- 
responding results without noise reduction are visually 
more noisy with identical growth exponents. 
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To proceed quantitatively we now introduce the dy- 
namic scaling ansatz which seems to describe well all 
our simulated results. We have studied the root mean 
square surface width or surface roughness (W), the av- 
erage mound size {R), the average mound height (H), 
and the average mound slope (M) as functions of growth 
time. We have also studied the various moments of dy- 
namical height-height correlation function, and these cor- 
relation function results (to be reported elsewhere) are 
consistent with the ones obtained from our study of W(t), 
R(t), H(t), and M(t). The dynamical scaling ansatz in 
the context of the evolving mound morphologies can be 
written as power laws in growth time (which is equivalent 
to power laws in the average film thickness) : W(t) ~ t@ 
; R(t) ~ t n ; H(t) ~ t K ; M(t) ~ t x ; £(i) - t^ z , where 
is the lateral correlation length (with z as the dy- 
namical exponent) and /3, n, K, A, z are various growth 
exponents which are not necessarily independent. In the 

steady state, when £(t) ~ L where L is the lateral sub- 
strate size, effective (3 vanishes as the surface roughness 
saturates to a value W S (L) = W(L, t — > oo) ~ L a , where 
a = (3z is the roughness exponent. We find in all our 
simulations n w z , and thus the coarsening exponent 
n, which describes how the individual mound sizes in- 
crease in time, is the same as the inverse dynamical ex- 
ponent in our model. We also find (3 — k in all our 
results, which is understandable in a mound-dominated 
morphology. In addition, all our results satisfy the ex- 
pected exponent identity (3 = K = n + A because the 
mound slope M ~ H / R. The evolving growth morphol- 
ogy is thus completely defined by two independent critical 
exponents (3 (the growth exponent) and n (the coarsening 
exponent), which is similar to the standard (i.e. without 
any diffusion bias) dynamic scaling situation where (3 and 
2 (= n~ l in the presence of diffusion bias) completely 
define the scaling properties. We note also that our neg- 
ative bias results (not shown here for lack of space) are 
completely consistent with the expected fl|,^2| Edwards - 
Wilkinson universality class 0j2l|,^] with our numerical 
findings being (3 = 0.25 and z = 2 in d=l+l, and (3 = 
(i.e. W ~ hit) and z — 2 in d=2+l. This is expected 
because our negative bias model explicitly introduces an 
inclination dependent downhill surface current [ p9| in the 
problem. 

In Fig. 3 we show our representative scaling results for 
nonequilibrium growth under surface diffusion bias con- 
ditions. It is clear that we consistently find ~ 0.5 in 
both d=l+l and 2+1 for growth under a surface diffusion 
bias. This (3 ~ 0.5 is, however, very different from the 
usual Poisson growth under pure random deposition with 
no relaxation where there are no lateral growth correla- 
tions. Our calculated asymptotic coarsening exponent n 
in both d=l+l and 2+1 is essentially zero (< 0.1) at long 
times. In all our results we find the effective coarsening 
exponent showing a crossover from n « 0.2 at early times 
to a rather small value (< 0.1) at long times - we believe 
the asymptotic n to be zero in our model. 
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FIG. 3. (a) The surface roughness W in 1+1 dimensions 
as a function of deposition time t in the L = 10000 system. 
Left inset: the growth exponent (3 calculated from the local 
derivative of log 10 W with respect to log 10 t. Right inset: Av- 
erage mound size as a function of deposition time, (b) The 
surface roughness in the 100 x 100 system. Left inset: the 
local growth exponent (3. Right inset: Average mound height 
vs time, (c) The average mound size in the 100 x 100 system. 
Left inset: the local coarsening exponent n calculated from 
the local derivative of log 10 R with respect to log 10 t. Right 
inset: the average mound slope vs time. 
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In comparing with the existing continuum growth 
equation results (of which there are quite a few) we 
find that none can quantitatively explain all our findings. 
Golubovic pH predicts (3 — 0.5, which is consistent with 
our simulations, but his finding of n = A = 1/4 in both 
d=l+l and 2+1 is inconsistent with our results (n < 0.1, 
A ~ 0.4 — 0.5) except in the limited transient regime. The 
analytic results of Rost and Krug [ljj also cannot explain 
our results, because they predict, in agreement with Gol- 
ubovic, that if (3 = 1/2, then n = A = 1/4, which is in 
disagreement with our findings. We also find our (3 to be 
essentially 0.5 independent of the actual value of Pjj—Pl, 
which disagrees with ref. 19. Our model obviously has no 
slope selection (A = 0), and therefore theories predicting 
slope selections do not apply. 

Before concluding we point out some reasons for why 
we believe our nonequilibrium limited mobility growth 
model to be a good zeroth order description for kinetic 
growth under a surface diffusion bias. The main rea- 
son for this is the success of the corresponding mini- 
mal growth model, introduced in ref. 23, in provid- 
ing a good zeroth order description of molecular beam 
epitaxial growth in the absence of any surface diffusion 
bias. The d=2+l critical exponents |^7| in the unbi- 
ased model || are (3 = 0.25 - 0.2 and a ~ 0.6 - 0.7, 
which are in quantitative agreement with a number of 
experimental measurements [ pl|j2^ ] where surface diffu- 
sion bias is thought to be dynamically unimportant. An 
equally significant theoretical reason is that the corre- 
sponding unbiased growth model |p3|-p7| is the only exist- 
ing nonequilibrium growth model which is known [ p2|j27| 
not to have the linear Edwards - Wilkinson 'V 2 h' term 
in its coarse-grained long wavelength continuum de- 
scription (an equivalent statement [^9| is that this is the 
only limited mobility model which has a vanishing surface 
current on a tilted substrate) - in fact, our negative bias 
model introduces precisely this £ V h' term by produc- 
ing a downhill current on a tilted substrate. (The other 
limited mobility nonequilibrium growth models [ j30|]3l| 
introduced in the literature are known to belong to the 
Edwards - Wilkinson universality class, and are therefore 
unsuitable for diffusion bias studies [p2|.) Therefore, the 
minimally biased version of this model which we study 
in this paper should be the appropriate zeroth order de- 
scription for nonequilibrium growth under surface diffu- 
sion bias conditions. Since an approximate continuum 
description for the unbiased model [^3| has recently been 
developed (26), one could use that as the starting point 
to construct a continuum growth model for the biased 
growth situation. Such a continuum description is, how- 
ever, extremely complex |2(| as it requires the existence 
of an infinite number of terms in the growth equation. It 
remains unclear that a meaningful continuum description 
for nonequilibrium growth under a surface diffusion bias 
is indeed possible even for our simple limited mobility 
growth model. 

We conclude by mentioning that any comparison to ex- 
periments, as has been done in several recent theoretical 



publications on this topic, is premature at this stage be- 
cause the real Ehrlich - Schwoebel barrier in experimen- 
tal systems Q is expected Q to be quite complicated, 
and simple growth models used by us and others fl-|2C|l 
most likely do not apply. What we have established in 
this paper is that even a very simple limited mobility 
nonequilibrium growth model leads to extremely complex 
dynamical interface morphologies under surface diffusion 
bias conditions. The fact that even a deceptively simple 
limited mobility growth model such as the one studied 
in this paper seems to defy a theoretical continuum de- 
scription is a strong indication that our understanding of 
nonequilibrium growth under a surface diffusion bias is 
far from complete. 

This work is supported by the US-ONR and the NSF- 
MRSEC. 
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